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Abstract 

The critical behavior at the frozen/active transition in the Domany-Kinzel 
stochastic cellular automaton (DKCA) is studied via a surface growth process 
in (1+1) dimensions. At criticality, this process presents a kinetic roughen- 
ing transition; we measure the critical exponents in simulations. Two update 
schemes are considered: in the symmetric scheme, the growth surfaces belong 
to the Directed Percolation (DP) universality class, except at one terminal 
point. At this point, the phase transition is discontinuous and the surfaces be- 
long to the Compact Directed Percolation universality class. The relabeling of 
space-time points in the nonsymmetric scheme alters the surface growth dra- 
matically. The critical behavior of rough surfaces at the nonchaotic/chaotic 
transition is also studied using the damage spreading technique; the exponents 
confirm DP values for the symmetric scheme. 
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I. INTRODUCTION 



The one- dimensional Domany-Kinzel stochastic cellular automaton (DKCA) is a com- 
pletely discrete system - temporally, spatially and in its state space - with applications in 
physics, chemistry, biology, computer science, etc. The DKCA also attracts interest 

as a particle system affording a test of certain conjectures regarding nonequilibrium critical 
phenomena ||. The DKCA has a unique absorbing ("vacuum") state; its phase diagram 
presents a critical line separating this absorbing phase from an active phase. Models with 
one absorbing state have been conjectured to belong generically to the directed percolation 
(DP) universality class [f|]. There is also good numerical evidence fU that the critical be- 
havior along the transition line in the DKCA belongs to the DP class, except at one of the 
terminal points, where the asymptotic behavior is known exactly and belongs to the compact 
directed percolation (CDP) universality class |)]||. At this terminal point the transition is 
discontinuous and we have in fact two absorbing states: the vacuum and the completely 
filled state. Using the damage spreading technique, Martins et al. JFj found another critical 
line separating the active phase into a nonchaotic and a chaotic phase. There is numerical 
evidence that the critical behavior along this transition line also belongs to the DP class, as 
expected on the basis of universality ||. 

The surface growth process generated by cellular automata (CA) was proposed by de 
Salles et al. [f| to study Wolfram's deterministic CA. These authors also used this process 
to identify the frozen/active phase transition in the DKCA [[[(J, where they showed that the 



Hurst exponent H attains a maximum at the transition. Recently, Atman and Moreira flT 



demonstrated that the growth exponent (3 W also presents a cusp at criticality, and is more 
appropriate for detecting phase transitions than the method of de Salles et al. They used it 
to construct the DKCA phase diagram and conjectured that the growth exponent method 
also can be used to detect phase transitions in other kinds of models. Recently, Redinz and 



Martins JT2| used the Hurst exponent method to find first and second order phase transitions 
in the g-state Potts model (for q — 1, 3, 5 and 10). Bhattacharyya [Oj studied the dynamic 
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critical properties of a related one- dimensional probabilistic cellular automaton, using two 
different procedures to generate the surface growth process. One of them is identical to the 
growth process studied by de Salles et al, and belongs to DP universality class. Thus, the 
conjectures presented by Bhattacharyya |TJ| which imply DP-like scaling for this model are 



also valid here. 

Basically, the procedure used by de Salles et al. |L0| and by Bhattacharyya ||13| , consists 



in transforming the spatiotemporal patterns generated by the DKCA to a solid-on-solid 
(SOS) particle deposition. At the critical line, it yields a surface growth process with kinetic 
roughening, and the critical exponents can be measured following the scaling concepts devel- 
oped by Family and Vicsek []14j . Very recently, similar methods were used by Lauritsen and 



Alava JT5]| to study the Edwards- Wilkinson equation with columnar noise, by Vespignani et 



al. to study sandpile models, and by Dickman and Munoz JT7J, to study the contact 



process (CP). 

In this work, we measure the scaling exponents at criticality in simulations, and com- 
pare them with known DP and CDP values. We use two different schemes to update the 
automaton, which lead to entirely different surface growth scaling properties. In Section II, 
we define the DKCA, describe the two update schemes, and show how the surface growth 
process is generated, at the frozen/active and nonchaotic/chaotic transitions. In Section 
III, we present our numerical results and discuss the values obtained, comparing them with 
the predictions for the scaling exponents proposed by Bhattacharyya Jl3| and Dickman and 
Munoz ]17|]. We discuss our conclusions in Section IV. 

II. DKCA SURFACE GROWTH PROCESS 
A. Model 

The DKCA was proposed by Domany and Kinzel who showed the existence of two 
phases: active and frozen. A more detailed study, using simulation, was performed by 
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Martins et al. 0, in which a new phase within the active region - a chaotic phase - was 
discovered through the damage spreading technique. 

The DKCA consists of a linear chain of L sites (i = 1, 2, L), with periodic boundaries, 
where each site i has two possible states, conveniently denoted by Uj = 0, 1. The state of 
the system at time t is given by the set {<7j(£)}. In contrast to the deterministic CA studied 
by Wolfram 0, the DKCA is probabilistic: the rules for updating the system are given by 
conditional probabilities, which depend on the neighbors. We study two different schemes, 
one symmetric, the other nonsymmetric. The symmetric scheme is the original one proposed 
by Domany-Kinzel Jl|, while the nonsymmetric was used Jl8|| to simplify the algorithm; its 



authors affirm that the two schemes are topologically equivalent. 

In the symmetric scheme, the process cxj(t) is defined on space-time points with i+t even. 
The state of site i at time t+1 depends on cjj_i(t) and via the transition probability 

P[<Ji(t+l)\<Ti-i(t),(Jitf_(t)], which takes the form: 
P(l|0,l) = P(l|l,0)= Pl , 

P(1|1,1)=P2 , 

P(1|0,0) = . 

Evidently, P(0|<7<_i, cr i+1 ) = 1 - P(l\a^i, a i+1 ). 

In the nonsymmetric scheme, the process is defined on all space-time points. The state 
of site i at time t + 1 depends on ovi(t) and Cj(t), rather than on o"j_i(t) and a i+ i(t), as 
in the symmetric scheme. The transition probability P[<7j(t + l)|o"j_i(i), o~i(t)} is identical to 
P [o"j (t + 1) | <7j_i (/!;), a^ilt)} given above for the symmetric case. 

It is easy to see that the two schemes are connected via a simple relabelling of space- 
time points (see Fig. 1). Consider, for example, a history h s in the symmetric scheme, 
that is, a sequence of configurations {erf (t)} for t = 0, ...,T (T finite), starting at t — 
with a finite number of active sites. Let the probability of this history, given the initial 
configuration, be P[/i 5 |{cr 5 (0)}]. A history h NS in the nonsymmetric scheme can be defined 
in the same manner. Note that there is a one-to-one correspondence between histories in 
the two schemes, given by 



(i) 



Since the transition probabilities in the two schemes are identical, the probabilities of 
corresponding histories are as well. To extend this correspondence to systems with pe- 
riodic boundaries, we note that if the nonsymmetric system has L sites, then the cor- 
responding symmetric one has 2L sites; in the mapping defined above, we now take 



An immediate result of this correspondence is that all scaling properties (e.g., critical 
exponents), as well as nonuniversal properties (e.g., phase boundaries between frozen, active, 
and chaotic phases in the pi — P2 plane), are identical in the two schemes. Corresponding 
histories naturally look different in the two schemes: the nonsymmetric scheme represents 
a rotating frame of reference in which, moreover, distances are rescaled by a factor of 1/2. 
Thus, for pi = 1/2 and P2 = 1, an interface between domains of l's and O's executes an 
unbiased random walk in the symmetric scheme, while in the nonsymmetric case such an 
interface has a mean velocity of 1/2. The "light cone" i — ±t in the symmetric scheme 
becomes the pair of lines i = and i — t in the nonsymmetric case. As will be seen below, 
this difference in frames of reference has important consequences for the surface dynamics 
in the nonsymmetric scheme. 

Depending on the values of the parameters (pi,P2), the asymptotic (t — > oo) state of 
the system is either frozen, with all sites having value 0, or has a finite fraction of sites 
with value 1, the active state. This is a second order phase transition, characterized by the 
critical exponents ofthe DP universality class. 



The surface growth process consists in accumulating (summing) all the values assumed 
by the variables <7j(r) over the first t time steps: 



i s = 2i NS -t (mod2L). 



B. Interface Representation 




(2) 



T = 
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The differences between the schemes become explicit at this point. In Figure 2 we show 
the temporal evolution of the automaton and the profiles generated by the accumulation 
method, close to criticality (p2 = 0.5, pi = 0.75), in each scheme. It is evident that the two 
schemes lead to entirely different profiles. (In this figure, we choose an initial condition of a 
single active site, to highlight the evolution of the automaton and profiles.) 

Thus, we obtain growth processes, the nature of whose correlations can be investigated 
through the analysis of the roughness w(L,t) fl9|], defined by 

w\L,t) = \[^{h % {t)-Kt)f) , (3) 

where h(t) is the mean value of hi(t) at time t, and the brackets < . . . > denote an average 
over realizations. 



We expect that w(L,t) has the scaling form 



w{L,t)~L°f(±;} , (4) 

where f(u) is a universal scaling function, a is the roughness exponent, z = ot//3 w is the 
dynamic exponent and j3 w is the growth exponent. The function f(u) = constant, at large 
times (t ^> L z ), and f(u) ~ u^ w at short times (t < So at short times, we expect 

w(t) ~ t@ w ; we measure j3 w from the slope of the log — log plot of w(L, t) versus t. At large 
times, the roughness saturates and becomes L-dependent: w(L, oo) ~ L a . The crossover 
time, £ x , between these two regimes grows as t x ~ L z . The exponents a and z are defined 
in the frozen phase just at the transition line. In the active phase, the roughness does not 
saturate, growing instead as w(L,t) ~ t 1 ^ 2 , corresponding to uncorrelated growth [[□]]. The 
relations above, used to measure the scaling exponents, are valid just at criticality. 

The profiles have self-affme properties quantified by the Hurst exponent H, defined via 

W(e) ~ e H , (5) 

where W(e) is the width of the interface on length scale e. We measure the Hurst exponent, 
H, in a profile generated very close to the transition. We apply the method introduced 



by Moreira et al f20|| , that consists in measuring the roughness around the straight line 
determined by a least-squares fit to a segment of the profile. The roughness W(L,e, t), at 
the scale e, is given by 

W(L,e,t) = ±5><(M) (6) 

i=l 

where the local roughness Wi(e,t) is defined by 

w*(e,t) = W) - +bj(e)]} 2 ■ 

j=i-e 

Oj(e) and bi(e) are the linear fitting parameters to the profile on the interval [i — e, % + e] 
centered at site i. 



C. Damage Spreading 

Martins et al. |7|] used the damage spreading technique to show that the active phase of 
the DKCA in fact consists of two phases, chaotic and nonchaotic. The order parameter of this 
transition is the difference between two replicas started with different initial configurations. 
One lets the system evolve until it attains a stationary state, and then a replica of the 
configuration is created with some sites altered (damage). The two replicas, one with state 
(Ti(t) and the other with state Q%{t), evolve with the same sequence of random numbers, and 
the difference between the configurations 

Ti(t) = \(Ti(t) -Qi(t)\ , 

is measured. The fraction of sites in the two replicas with <7j 7^ Qi is called the Hamming 
distance, defined as 

The stationary Hamming distance is null in the non-chaotic phase and positive in the chaotic 
phase. 



S 



To study the chaotic/non-chaotic boundary, we use a slightly different method, where 
the difference between the two automata is used to generate the surface growth process, as 
we did in the accumulation method 

M*) = Er<(r). (7) 

r=0 

Thus, the profile generated by the difference between the replicas behaves exactly as the 
profiles generated in the frozen/active boundary: the roughness reaches a stationary value 
in the non-chaotic phase and grows indefinitely in the chaotic phase. This behavior can be 
understood if we note that the difference between the replicas vanishes in the non-chaotic 
phase, which implies no contribution to the height hi(t), and is positive in chaotic phase, 
implying steady growth in the height. 

In order to preserve the stationary density of active sites we generate a "rotation" damage 
at a certain time to, in which the replica is rotated 180° with respect to the original system, 
that is g(i,to) = cr(i + L/2,t ), subject to the periodic boundary condition. 



D. Theoretical Descriptions 

Theoretical descriptions of surface growth scaling at absorbing-state phase transitions 
were proposed by Bhattacharyya and Dickman and Munoz fl7fl . Bhattacharrya pro- 



posed an analytical treatment in analogy to the random deposition (RD) process: considering 
an initial disordered state, the growth process can be described by a continuum equation, 
very similar to the RD process: 

^^- = F + r,(x,t) [RD] (8) 

where F is the average number of deposited particles and r)(x,t) corresponds to the white 
noise (< r](x,t) >= 0) in the deposition. 

The difference between the surface growth processes generated by DKCA and RD lies in 
the noise correlations. While RD involves spatially and temporally uncorrelated noise, the 
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correlations in time and space developing in the DKCA appears in the noise fluctuations of 
the accumulation method. For values of (pl,p2) away from the critical line (in the active 
phase), the correlation length £ and correlation time r of the DKCA are finite, which means 
that the noise in the deposition process is correlated over short ranges. In this limit the 
noise autocorrelation decays exponentially [|nj: 

< V (x, t),v(x', t') >~ e -W-x'\/t e -\t-t'\/T _ (9) 
Thus, the noise appears uncorrelated for times greater than r, and the RD exponents 



are obtained in this limit. This behavior was confirmed in earlier simulations ||11|| . As we 
approach the transition line, £ and r increase, and it takes longer for the growth process to 
reach the RD limit. Finally, at the critical line, both £ and r diverge, and the correlations 



are long ranged, represented by a power law decay of the noise autocorrelation |19 



< r)(x,t),r)(x>,t') >~ \x-x'\- 2(3/ ^\t-t'\- 2f3/ ^ , (10) 

where 3, v± and u\\ are, respectively, the critical exponents for the order parameter, corre- 
lation length and correlation time of the DKCA. 

The value of the growth exponent 3 W at the critical line can be derived from the con- 
tinuum equation (f|) and the noise (|i~0D ; the width for an infinite substrate is expected to 



increase as the following power of the time |L3 



w{oo,t) -t 1 -^ . (11) 
As it is believed that this transition belongs to the DP universality class, the value of growth 



exponent is expected to be [21 



3 0.27649 

8 W = 1- — = 1 ~ 0.8405 . 

h v\\ 1.733825 

In previous work |IT| , Atman and Moreira showed that 3 W attains a maximum at the 
phase transition (see Fig.3), and measured its value along the transition line of the DKCA. 
This behavior of the exponent 3 W in the vicinity of the phase transition can be understood as 
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follows. The growth rate, dhi/dt, at site i is proportional (in the frame of reference moving 
with the average velocity d < h > /dt) to the excess activity at that site. Away from the 
critical point, the activity has a finite correlation lenght £ and correlation time r. Thus on 
scales much greater than £, r, the noise driving the surface growth is uncorrelated, and this 
process falls in the RD class, with (3 W = 1/2. At the critical point, by contrast, £ and r 



diverge and we have instead the scaling relation (3 W = 1 — 6 \ T% , where the exponent 6 is 
defined through the relation p(t) ~ t~ e for the initial decay of the activity density p, at the 
critical point, starting from p(0) = 1. Since, in one dimension, 1 — 6 > 1/2, we expect a 
jump in (3 W at the phase boundary. In simulations of finite-sized systems, we expect not a 
discontinuity in (3 W but a sharp peak at the transition - see Figure 3 (very near to the critical 
point, £ > L, so that independently fluctuating regions are not present in the simulation. 
Below the transition, the apparent value of f3 w — > due to the short lifetime of the activity). 
It is interesting note that for the DP universality class, 1 — 9 ~ 0.55, 0.27 and for d = 2, 3 
and 4, respectively. Thus we should expect (3 W to decrease at the phase boundary in d = 3. 

Since the saturation of surface width is forced by the DKCA, the crossover time t x 
behaves exactly as in DP, and the dynamic exponent is given by: 

t x ~ L z , z = z DP = v\\/v± ~ 1.5808 . (12) 

Thus, the roughness exponent at the criticality is given by 

a = z ■ p w ~ 1.3286 . (13) 



Dickman and Munoz [17 studied the contact process using the surface growth repre- 



sentation. They demonstrated that the Hurst exponent shows clear signs of anomalous 



scaling (a > H) , but no evidence of multiscaling. They verified Lopez' scaling relation |22" 
H = a — zk. Here, k is the exponent associated with the divergence of the mean-square 
height gradient in the continuum growth equation that describes the contact process and 
related models, such as DP. Inserting the known values in the relation above we have for 
the Hurst exponent in the DP class: 
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H = a — zk ~ 0.643 



(14) 



III. RESULTS 

In Table 1, we summarize our results for the scaling exponents at the frozen/active and 
nonchaotic/chaotic transitions, and compare them with the values for the DP and CPD 
universality classes. 

To extract the exponent values from our simulation data, we used the relations w(L, oo) ~ 
L a , valid at large times, w(L, t) ~ t^ w , valid at short times and t x (L) ~ L z . The results show 
a strong dependence on the scheme used - symmetric or nonsymmetric. In the simulations, 
we average 10000, 5000, 2500, 1000, 500, 250 and 100 samples at the critical point {p c 2 ,p\) 
in systems with L = 50, 100, 200, 500, 1000, 2000 and 5000 sites, respectively. The initial 
condition in these samples was random, with 50% of sites active. 

The Hurst exponent was measured following the procedure explained in subsection 2.2. 
The results, shown in Table I, represent an average over 100 random initial configurations in a 
system with L = 10000. We observe a significant change in the value of the Hurst exponent 
depending on the scheme used to update the automaton: in the symmetric scheme, we 
have H > 1/2, denoting a positive correlation in the profile; in the nonsymmetric scheme 
H = 0.25(3), denoting a negative correlation. This behavior can be understood considering 
the nonsymmetric scheme as a deposition over a moving reference frame, which implies a 
lateral propagation of correlations. 

A. Symmetric scheme 

For the symmetric scheme, our results for the critical exponents agree with the DP values, 
except at terminal point p 2 — 1 where CDP values were obtained. The critical points (p^Pi) 
were determined through the growth exponent method |Tl| (Figure 3). This method consists 
in fixing p 2 and varying pi until the maximum of the growth exponent f3 w is attained. 
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To determine the crossover time t x (L), we plot the fraction of realizations with at least 
one active site as a function of time (see Fig. 4), and define the crossover time such that 
half of the initial sample has frozen. The inset of Figure 4 shows the power law behavior of 
the crossover time; the slope of this line corresponds to the exponent z. 

To obtain the saturation width, we let all samples evolve to the absorbing state, at a given 
system size, and determine the final averaged roughness. The exponent a was measured as 
the slope of w(L, oo) versus L in a log-log plot. In Figure 5 we present the results for the 
saturation roughness in the cases: directed percolation (p 2 = 0.5, pi = 0.749, symmetric), 
directed percolation (p 2 — 0,pi — 0.8095, damage spreading), compact directed percolation 
(p2 = 1, Pi = 0.5, symmetric) and nonsymmetric DP (p 2 = 0.5, p\ = 0.749, nonsymmetric). 

In order to verify Family- Vicsek scaling, we use the measured values for the scaling 
exponents to collapse the width curves at different system sizes to a single curve, as shown 
in Figure 6. Note the collapse of the width curves, corroborating the Family- Vicek scaling 
relation. 

It is important to note that the exponents f3 w and z measured for the chaotic/nonchaotic 
transition at p 2 = are slightly different from the exponents measured away from this 
point. This represents evidence of long-range correlations due the coincidence of the damage 
spreading and frozen/active transitions at this point, as pointed out by Grassberger ||. 

B. Nonsymmetric scheme 

A significant change in the roughening occurs in the nonsymmetric scheme, as shown in 
Figure 7; we observe two distinct regimes in the roughness growth: a strongly correlated 
regime, for times t ~ L, and a weak correlation regime, for longer times. Again, this 
behavior can be understood by considering the nonsymmetric scheme as deposition in a 
moving reference frame. The correlations inherent in the dynamics are propagated by the 
moving reference frame until they reach the system size; then the correlations due the local 
rules of the automaton take over, decreasing the growth rate of the roughness. These two 
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regimes for the roughness growth implies that is not possible colapse all the curves using 
the Family- Vicsek scaling law. 

The apparent exponent values measured in the nonsymmetric scheme are markedly 
smaller: a ~ 0.93(2) for nonchaotic/chaotic and frozen/active (at p 2 ^ 1) transitions; 
at p 2 = l,pi = 0.5, a ~ 0.984(7). As discussed above, the growth exponent presents two 
values, depending on the roughness growth regime. The dynamic exponent z must assume 
the DP value, as discussed in Sec. 2.1; in fact, z ~ 1-6(1) for nonchaotic/chaotic and 
frozen/active (at p 2 ^ 1) transitions, while z ~ 1-9(1) at p 2 = l,pi = 0.5. 

IV. CONCLUSIONS 

Growth surfaces generated by the spatiotemporal patterns of the DKCA along its critical 
lines are studied. The critical roughening exponents, expected to belong to the DP univer- 
sality class, were measured using power law relations valid at criticality. Except for the 
terminal point p 2 = 1, all the scaling exponents agree with the DP values, in the symmetric 
scheme, and the scaling law j3 w = a/z remains valid. At p 2 = 1, we confirm CDP values 
for the exponents. Since the fluctuations in uncorrelated regions are effectively superposed, 
it is not surprising that the apparent values of f3 w and a are smaller in the nonsymmetric 
scheme. At the nonchaotic/chaotic transition, the exponents measured also agree with the 
DP values. 

Acknowledgements We thank the Brazilian agencies CNPq and Fapemig for financial 
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TABLES 



TABLE I. Summary of the scaling 


; exponents values, 


d=l. 




Previous work 


a 


P w 


z 


H 


DP 


1.3286 


0.8405 


1.5808 


0.643 


CP (Simulational) |7) 


1.33 


0.839(1) 




0.63(3) 


CA (Simulational) @ 




0.837(11) 






CDP 


2 


1 


2 


1 


Present work - Symmetric scheme 


frozen/active P2 = 0.5 


1.32(1) 


0.82(2) 


1.59(1) 


0.61(3) 


frozen/active P2 = 1 


2.01(1) 


0.99(1) 


2.08(5) 


0.99(2) 


nonchaotic/chaotic p\ = 1 


1.325(9) 


0.81(1) 


1.61(1) 


0.60(3) 


nonchaotic/chaotic p2 = 


1.32(1) 


0.78(2) 


1.64(2) 


0.61(3) 
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Figure Captions 

Figure 1 Spatial representation of DKCA, in symmetric (left) and nonsymmetric (right) 
schemes, showing that the spatio-temporal patterns are identical in the two schemes, i.e., 
corresponding histories are identical. 

Figure 2 Growth surfaces generated by spatiotemporal patterns of the DKCA with different 
update schemes: the symmetric scheme is shown on the left and the nonsymmetric on the 
right. Upper panels: temporal evolution of DKCA. Black sites are active; time increases 
upward. Lower: profiles generated by the accumulation method; the fill color is changed 
every 50 steps. System size is L = 500 for the symmetric and L = 250 for the nonsymmetric; 
900 time steps are shown. Both systems are very close to criticality (p 2 = 0.5, p 1 = 0.75) in 
the active phase. 

Figure 3 Growth exponent f3 w in the DKCA for several system sizes in the symmetric 
scheme. Note that f3 w attains a maximum at the frozen/active transition and depends 
strongly on the system size. The transition point was chosen as the f3 w value at system size 
L = 10000, where we observe a sharp transition. In this example, p 2 = 0.5 and p\ = 0.749. 
Figure 4 Dynamic exponent z for the DKCA interface representation. The density p(t) 
of active samples in function of time for several system sizes is shown, for the symmetric 
scheme (p 2 = 0.5, pi = 0.749). The horizontal line highlight the value p = 1/2, which 
corresponds to the crossover time. The inset shows the crossover time t x in function of the 
system sizes. The slope of this curve is the value of the dynamical exponent z. The error 
bars are calculated considering an error of 1% in the number of samples froze at a given 
time. 

Figure 5 Roughness exponent a for the DKCA interface representation. Four cases are 
shown: directed percolation (p 2 = 0.5, p\ = 0.749 symmetric), directed percolation (p 2 = 0, 
Pi = 0.8095 damage spreading), compact directed percolation (p 2 — 1, p 1 — 0.5 symmetric) 
and nonsymmetric DP {p 2 = 0.5, p\ = 0.749 nonsymmetric). The line is the power law 
regression for the data and furnishes the value of the roughness exponent a. The error bars 
are the standard deviations of the saturation width over realizations at each system size. 
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Figure 6 Family- Vicsek scaling. Upper panel: the width of the generated profiles at different 
system sizes. Lower: collapse of the curves above using the exponent values measured 
through numerical simulation at the criticality (j>2 = 0.5, pi = 0.749). 

Figure 7 Profile roughness behavior for nonsymmetric scheme. Note the two regimes in 
the w(L, t) x t curve; for t ~ L we have a srongly correlated regime with (3 W ~ 0.65 and for 
L < t < t x , a weakly correlated regime with f3 w ~ 0.45. 



19 



0.80 



0.70 



0.60 



0.50 



0.40 



0.70 



p 1 =0.749(1) [transition] p 2 =0.5 




L=10000 

L=5000 

L=2000 

L=1000 

L=500 

L=200 

L=100 



0.75 



0.80 
Pi 



0.85 



0.90 



10 5 



10' 



p 2 =0.5 
p 1= 0.749 




10 



10 u 



10' 



1(T 



10 c 



m 

CO 



10 



10 




10 



10 



10 



t/L 



1.567 



10 



e L=100 
a L=200 



o o L=500 

a a L=1000 

o < L=2000 
v v L=5000 



> L= 10000 



10' 



; 


1 ; 

weakly 




correlated 




regime , ^ ~_ 


: / 

: strongly / 


/ / - 

/ / " : 
/ ' ; 


correlated / x 




regime / X 


nonsymmetric scheme 




p 2 =0, p^O.8095 




L=50 


/ x 
/ / 
/ / 


L=500 




L=5000 : 



10° 10 2 10 4 10 6 10 



